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Abstract 

We use variable-pressure neutron and X-ray diffraction measurements to determine the uniaxial and bulk compressibilities of 
nickel(II) cyanide, Ni(CN) 2 . Whereas other layered molecular framework materials are known to exhibit negative area compress¬ 
ibility, we find that Ni(CN )2 does not. We attribute this difference to the existence of low-energy in-plane tilt modes that provide 
a pressure-activated mechanism for layer contraction. The experimental bulk modulus we measure is about four times lower than 
that reported elsewhere on the basis of density functional theory methods [Phys. Rev. B 83, 024301 (2011)]. 
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1. Introduction 

From a lattice dynamical perspective, molecular frameworks 
have provided an extraordinary number of interesting exam¬ 
ples that challenge our understanding of how materials should 
respond to external stimuli such as temperature and pressure. 
Negative thermal expansion in the Zn(CN )2 family d |2l, 
pressure-induced and thermal amorphisation in metal-organic 
frameworks (MOFs) such as ZIF -8 Eii, and “breathing” tran¬ 
sitions in the MIL-53 and MIL-101 systems are all popular ex¬ 
amples EEl. A related phenomenon to receive particular re¬ 
cent interest is that of negative compressibility, whereby a ma¬ 
terial actually expands in one or more directions under appli¬ 
cation of hydrostatic pressure dUl. While it is now emerging 
that negative linear compressibility (NLC) might be somewhat 
more widespread amongst molecular frameworks and molec¬ 
ular crystals than originally thought 0 , the much rarer prop¬ 
erty of negative area compressibility (NAC)—where a system 
reduces its volume by expanding simultaneously in two orthog¬ 
onal directions—remains known to occur in only a handful of 
materials |[T0l[TTl[T2l. The development of a microscopic un¬ 
derstanding of NAC and its lattice dynamical origins is crucial 
if we are ever to design improved NAC candidates for applica¬ 
tion in next-generation ferroelectric sensor technology Q . 

Arguably the clearest lattice dynamical signature of NAC es¬ 
tablished to date is the existence of so-called Lifshitz modes 
in layered framework materials. These modes are low-energy 
vibrations of weakly-coupled layers in which atomic displace¬ 
ments are polarised perpendicular to the layer plane [Fig. 

The relationship with NAC arises from a combination of a large 
and negative elastic compliance S 13 —linking interlayer com¬ 
pression with layer expansion—and a small value of the elastic 
constant C 33 which refiects facile compression along the layer 
stacking axis. As a consequence, one design strategy for identi¬ 



Figure 1: The Lifshitz instability in layered materials involves displacement 
of atoms within layers in a direction parallel to the layer stacking axis. In¬ 
creased population of Lifshitz modes results in a reduction of the effective 
cross-sectional area of each layer, which couples to an expansion of the inter¬ 
layer separation. The Lifshitz mode is associated with negative area compress¬ 
ibility in cases where volume reduction under hydrostatic pressure results in 
rapid inter-layer collapse that couples with damping of the Lifshitz modes. 

fying new NAC candidates has been to focus on layered systems 
with very weak inter-layer interactions, such as is observed in 
the key NAC system silver(I) tricyanomethanide ifTTl . 

It was in this context that we sought to study the high- 
pressure structural behaviour of nickel(II) cyanide, Ni(CN) 2 . 
The structure of this material is described by stacked square 
layers, each assembled from square-planar Ni^^ ions connected 
via linear Ni-C-N-Ni linkages [Fig. |[T3l [141. The inter¬ 
action between layers is considered to be weak on the basis 
of a number of observations. First, Ni(CN )2 exhibits strong 
positive thermal expansion (PTE) behaviour along the layer 
stacking axis that is coupled to area negative thermal expansion 
(NTE) within the layer plane IH; note the magnitude of ther¬ 
mal expansion is usually inversely proportional to the strength 
of bonding interactions ca Second, there is no long-range pe¬ 
riodicity within the stacking arrangement, a structural feature 
which has also been associated with weak inter-layer interac¬ 
tions ca. And, third, the phonon dispersion curves determined 
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Figure 2: Representation of the structure of Ni(CN) 2 . (a) Square-grid layers 
are stacked along the c crystal axes. Successive layers are shifted by |<100> 
lattice vectors. The existence of choice in this shift direction results in a loss of 
long-range periodicity along the stacking axis, (b) A single layer is assembled 
from square-planar Ni^^ cations (shown here in teal) connected via linear Ni- 
C-N-Ni linkages (C/N atoms shown in grey). The orientation of the cyanide 
ions is disordered throughout the structure. The view direction in (b) is parallel 
to the c axis. 

using density functional theory (DFT) calculations reflect the 
existence of Lifshitz-like excitations ca. So the established 
design criteria for NAC in molecular frameworks would suggest 
that Ni(CN )2 should be a prime candidate in which to observe 
the unusual phenomenon of NAC. 

Here we apply a combination of variable-pressure neutron 
and X-ray powder diffraction measurements to determine ex¬ 
perimentally the compressibility behaviour of Ni(CN )2 over the 
pressure range 0 < p < 2.5 GPa. We And no evidence for 
NAC behaviour, a result we rationalise in terms of the existence 
of low-energy modes with negative Griineisen parameters that 
compete with the Lifshitz modes in allowing volume collapse 
via layer densiflcation. Subsequent to carrying out our mea¬ 
surements, a related study has appeared in the preprint literature 
(m. The experimental results reported by us and by Ref.fTTlare 
broadly similar; however our analysis differs in some important 
respects. In particular, we interpret our compressibility data in 
the context of the Griineisen parameters associated with differ¬ 
ent types of NTE modes supported by the Ni(CN )2 structure. 
Our own paper is arranged as follows. We begin by describing 
the experimental methods used in our study. We then present 
the pressure-dependent lattice parameter variation we observe 
experimentally and report the corresponding lattice and bulk 
compressibilities. Our manuscript concludes with a discussion 
of possible microscopic mechanisms responsible for this be¬ 
haviour, placing our results in the context of previous studies 
of Ni(CN )2 and of other NAC materials. 

2. Experimental Methods 

2.7. Synthesis 

Powder samples of Ni(CN )2 were prepared according to 
one of two methods. The first method involved heating 
Ni(CN)24H20 (Alfa Aesar, >99%) under vacuum at 200 °C for 
24 h. The dehydration process is coupled to a change in colour; 
the yellow product was stored in a vacuum desiccator to prevent 
rehydration. The second method involved reaction of aque¬ 
ous solutions of nickel(II) nitrate (Sigma Aldrich, >99%) and 
potassium tetracyanonickelate(II) (Sigma Aldrich, >99%). The 


precipitate formed was isolated by vacuum filtration, washed 
(H 2 O) and dried to afford Ni(CN )2 • 4 H 2 O. This product was 
dehydrated and stored as described above. 

2.2. Variable-pres sure neutron diffraction 

Neutron diffraction measurements were carried out using 
the PEARL instrument at the ISIS neutron spallation source. 
Hydrostatic pressure was applied using a V4 Paris-Edinburgh 
press. A powdered sample of Ni(CN) 2 , prepared as described 
above, was loaded into a TiZr gasket together with a Pb pel¬ 
let and fluorinert as pressure-transmitting medium. The gas¬ 
ket was positioned within a zirconia-toughened alumina anvil 
with a Los Alamos single toroidal profile. The sample was sub¬ 
jected to a series of eight force values that corresponded to hy¬ 
drostatic pressures 0 < p < 0.78 GPa. Pressures were calcu¬ 
lated using the third-order Birch-Mumaghan equation of state 
known for Pb (Vq = 30.354A^ Bq = 42.0GPa, B' = 5.71). 
Neutron scattering data were collected over the J-spacing range 
0 < J < 4.1 A. Constrained Rietveld refinements were used to 
extract unit cell parameters via the GSAS refinement package 
ifT^ . A crystallographic approximant to the “disorder stack” 
model reported in Ref. [14] was used as the basis of these re¬ 
finements (I4i/amd,a ^ 4.9 A, c ^ 12.8 A), together with a 
stacking-fault model as implemented in GSAS ifTSll . Our data 
were not sufficiently discriminating to allow robust refinement 
of positional parameters. 

2.3. Variable-pres sure synchrotron X-ray diffraction 

Synchrotron X-ray diffraction measurements were carried 
out using the 115 beamline at the Diamond Light Source. A 
powdered sample of Ni(CN) 2 , prepared as described above, was 
loaded under a controlled atmosphere into a modified Merrill- 
Bassett diamond-anvil cell with diamond culets of 450 pm. Sil¬ 
icone oil was used as pressure-transmitting medium. Ruby 
chips, evenly distributed in the pressure chamber, were used to 
measure the applied pressure according to the ruby fluorescence 
method |[T9j[20l. The variation in pressure values determined 
from different chips was found to be less than 0.1 GPa. For our 
measurements, five pressures over the range 0 < /? < 2.5 GPa 
were applied. The observed intensities were integrated as a 
function of 26 using the software FIT2D 1211 [22]| in order to 
give one-dimensional diffraction profiles. A powder sample of 
Ce 02 was used to calibrate the beam centre and to determine 
the sample-to-detector distance. As for the neutron diffraction 
data, constrained Rietveld refinements were used to extract unit 
cell parameters via the GSAS refinement package im. The 
same structural model described above was used as the basis of 
these refinements, and again our data were not sufficiently dis¬ 
criminating to allow robust refinement of positional parameters. 

3. Results 

The time-of-flight neutron diffraction patterns and corre¬ 
sponding fits obtained from Rietveld refinement are shown in 
Fig. The coexistence of contrasting sharp and broad reflec¬ 
tions noted elsewhere C3 as evidence for stacking disorder 
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Figure 3 : Variable-temperature time-of-flight neutron diffraction patterns col¬ 
lected for Ni(CN)2 together with fits obtained from Rietveld refinement. Data 
are shown as black markers; fit as a red line, and difference (data - fit) in blue. 
Tick marks are given for the various phases included in the fit: (bottom-top) 
Ni(CN)2, Pb, AI2O3 and Zr02 components of the zirconia-toughened alumina 
(ZTA) anvils, and (optionally) the high-pressure Ni(CN)2 phase. 


is clearly visible in our newly-collected data and appears im¬ 
mune to the application of hydrostatic pressure. In contrast to 
the finding of Ref. [T7] we do not find any evidence of a struc¬ 
tural phase transition below 0.1 GPa: all data sets collected over 
the pressure range 0 < ;? < 0.5 GPa are well fit using a struc¬ 
tural model based on the ambient-pressure structure d. At 
pressures above 0.5 GPa we find a set of additional refiections 
that can be indexed (very tentatively) on the basis of a new 
Ni(CN )2 cell consistent with a defect Prussian blue architecture 
(Im3m, a ^ 4.8 A). In the absence of any noticeable disconti¬ 
nuity in the compressibility of the ambient phase we consider 
these additional peaks to arise from a slow first-order transition 
to a high-pressure phase that may be the p > 0.7 GPa phase pos¬ 
tulated in Ref. [TtI Given that only a handful of relatively-weak 
peaks contribute to this minor component of the diffraction pat¬ 
tern, we cannot have any real confidence in the structural model 
we have used for this phase in our Rietveld refinement. Nev¬ 
ertheless, the (meta-)stability of the ambient phase to higher 
pressures allows us to continue to study its compressibility be¬ 
haviour over the entire pressure regime 0 < p <2.5 GPa. 

The variation in lattice parameters extracted from our Ri¬ 
etveld fits to these neutron diffraction data—together with 
those for the X-ray diffraction data (themselves not shown)—is 
shown in Fig.[^ What is immediately obvious is that the struc¬ 
tural response to hydrostatic pressure is dominated by rapid col¬ 
lapse of the inter-layer separation, which corresponds to the 
c lattice parameter. By contrast the a parameter, which de¬ 
scribes the Ni.. .Ni separation within an individual square sheet. 



;j(GPa) 


Figure 4 : Compressibility behaviour of Ni(CN)2. (a) Relative change in lattice 
parameters under hydrostatic pressure: the mechanical response of the system 
is dominated by compression along c. (b) The corresponding linear compress¬ 
ibilities extracted from the empirical fit to the data in (a) according to Eq. 

(c) The volume equation of state, shown with both second-order {B' = 4 ) and 
third-order Birch-Murnaghan fits as described in the text. 


changes remarkably little under hydrostatic pressure. Of partic¬ 
ular note is the absence of any increase in this parameter with 
pressure, which rules out NAG behaviour in Ni(CN) 2 . 

The lattice parameter data were fitted in the program PAscal 
||23]| according to a set of empirical expressions of the form 

€{p) = (Q + a{p- (1) 

as described in Ref. [24| The corresponding compressibilities, 
obtained via the isothermal pressure derivative of Eq. Q, are 
shown in Fig. [Jb). Extreme deviations from linear behaviour 
are observed for the c axis compressibility at low pressures, 
which suggests strongly anharmonic inter-layer interactions. 
The linear compressibilities extracted from these fits over the 
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entire pressure range 0 < p <2.5 GPa are 

Ka = +1.8(2) TPa“\ (2) 

= +45(2) TPa“'. (3) 

Hence Ni(CN )2 is essentially incompressible within the layer 
directions and is about as strongly compressible in a direction 
parallel to the layer stacking axis as the NAC material silver(I) 
tricyanomethanide {K^ = +66(20) TPa“^) ifTTI . 

A third-order Birch-Murnaghan equation of state Ga de- 
scribes well the pressure-dependence of the unit cell volume 
[Fig. I^c)]. The corresponding zero-pressure bulk modulus is 
determined as = 11.7(1.8) GPa with a first pressure deriva¬ 
tive B' = 12(3). Taken together these values refiect a compli¬ 
ant material which rapidly stiffens on application of hydrostatic 
pressure 1^ . Our values are consistent with other cyanides and 
layered molecular frameworks ifTTl l27l but contrast with those 
reported in Refs. dSl [HI = 105(2) GPa (0 < p < 0.1 GPa; 
experiment), 15.4(2) GPa (O.l < p < 18 GPa; experiment), and 
63.4 GPa (DFT). 

4. Discussion and Conclusions 

There are two obvious questions posed by these results, and 
we proceed to answer these in turn. The first is why there might 
be such a discrepancy between the compressibility values de¬ 
termined here and those reported in Refs. [T6| and [T^ Whereas 
the second—and arguably the more important—concerns why 
Ni(CN )2 does not show NAC when its structure possesses all 
the design motifs thought to favour the phenomenon. Answer¬ 
ing this particular question is key to informing future design 
strategies for targeting NAC behaviour in other materials. 

With respect to the apparent discrepancies in compressibility 
amongst the various studies of the lattice dynamics of Ni(CN) 2 , 
we address first the difference between experimental and com¬ 
putational (DFT) results. The DFT study of Ref. |T6| was al¬ 
ways going to be difficult because the layer stacking sequence 
in Ni(CN )2 has no true periodicity along the stacking axis d, 
and so the material structure is inconsistent with the presump¬ 
tion of crystallinity inherent to DFT methods. The work-around 
devised in that study was to use a modified layer stacking ar¬ 
rangement in which the layers were separated by twice their 
normal distance and placed directly above one another rather 
than offset by ^(100) as known to occur in the real material. 
This meant that cell relaxation was not possible and the all- 
important compressibility in the c direction determined compu¬ 
tationally is unlikely to have any physical meaning. Even if a 
different stacking arrangement had been implemented, it would 
likely have proven difficult for DFT to determine the correct 
compressibility given the non-negligible role of thermal fluc¬ 
tuations in modifying the elastic properties of cyanide frame¬ 
works with (strictly) linear M-C-N-M linkages 1^ . 

The experimental data of Ref. [17] are much more consistent 
with our own results, with the exception of the compressibility 
behaviour at very lowest pressures 0 < p < 0.1 GPa. It was 
reported that, over this pressure regime, very little change in 


Table 1: Comparison of compressibility behaviour of Ni(CN )2 determined us¬ 
ing neutron and X-ray diffraction methods (our study) and neutron diffraction 
methods alone (Ref. E). 

This study Ref.fTTl 


Ka (TPa-‘) 

+1.8(2) 

+2 

+3.5 

Kc (TPa-i) 

+45(2) 

+3.5 

+45 

Bo (GPa) 

11.7(1.8) 

105(2) 

15.4(2) 

B' 

12(3) 

= 4 

= 4 

Pressure range (GPa) 

0-2.5 

0-0.1 

0.2-2.0 


the unit cell dimensions is observed but application of hydro¬ 
static pressure (Ar gas) induces an irreversible shift in cyanide 
stretching frequencies that the authors of Ref. [T7] associate with 
a transition in cyanide ordering. This process is difficult to ra¬ 
tionalise in terms of the extremely small pAV driving force (ca 
20 Jmor^) and the absence of any variation in neutron scatter¬ 
ing intensities as would be expected with the C/N sensitivity of¬ 
fered by the technique. We also find no evidence in our neutron 
scattering patterns for a difference in cyanide order between 
the ambient and high-pressure structures of Ni(CN) 2 . Given 
that the neutron and Raman scattering measurements of Ref.fTTl 
were performed under different experimental conditions, and 
given the absence of a pressure marker in the low-pressure neu¬ 
tron scattering study, we suggest there is sufficient ambiguity 
regarding this low-pressure regime that the longer-range com¬ 
pressibilities over the pressure range 0 < p < 2.5 GPa more 
likely refiect the true pressure response of Ni(CN )2 [Table 
Having established that Ni(CN )2 shows conventional (posi¬ 
tive) area compressibility in the plane perpendicular to its stack¬ 
ing axis, we proceed to discuss how this observation might be 
consistent with the existence of Lifshitz modes in its lattice dy¬ 
namics. These modes are characterised by a positive Griineisen 
parameter 


V (doj\ 


(4) 


where oj is the mode frequency and V the unit cell volume Q 
Consequently as V decreases with increasing hydrostatic pres¬ 
sure, the value of oj increases and so the mode occupation num¬ 
ber n(oj) ^ k^TIhco decreases; this change in frequency is ob¬ 
vious in the spectroscopic measurements of Ref. [TT] In other 
words, the degree of layer rippling is reduced at high pressures 
and the layer cross-sectional area A = might be expected to 
expand as a result. 

But the Ni(CN )2 structure supports an additional type of low- 
energy phonon mode in which [Ni(C/N) 4 ] units rotate about 
an axis parallel to the stacking axis [Fig. g USED. These 
in-plane tilts give rise to C/N displacements that are polarised 
within the layer plane and are related to the NTE mechanism 


^Note that the sign of y is misassigned in Ref. [16] (likely as a result of the 
particular stacking approximation used), which explains why that study predicts 
volume NTE at low temperatures in contrast to the PTE behaviour observed 
experimentally na. 
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Figure 5: Representative in-plane tilt mode responsible for both area-NTE and 
positive area compressibility. 

thought to operate in the Prussian Blues (301 [311. The corre¬ 
sponding Gnineisen parameters are now negative because pop¬ 
ulation of these modes leads to volume reduction: while the 
a lattice parameter decreases as for the Lifshitz modes, this 
effect is no longer counterbalanced by a change in the stack¬ 
ing axis repeat c. Because y < 0 the energy of these modes 
will decrease on application of pressure— i.e. their occupation 
numbers n{oj) will increase—leading to a contraction in cross- 
sectional area A and hence positive compressibility within the 
(a, b) plane. The absence of NAC in Ni(CN )2 would suggest 
that this second mechanism outweighs the Lifshitz mechanism 
for this material, but what is really needed for quantitative in¬ 
sight here is a computational study that takes into account the 
experimentally-observed stacking arrangement and allows in¬ 
terpretation of the inelastic neutron scattering results presented 
elsewhere lEl. One possible avenue would be the development 
of an ab initio force field that might then be used to drive finite- 
temperature molecular dynamics simulations, as has recently 
been carried out for Zn(CN )2 (321 • 

In terms of materials design strategies for engineering new 
NAC systems, what are the lessons suggested by our study? 
Perhaps the most obvious is that the existence of Lifshitz modes 
in a layered material does not guarantee NAC. We have sug¬ 
gested elsewhere that a natural domain in which to search 
for NAC candidates is amongst materials that show area-NTE 
Guia. What has emerged in our study is that the existence of 
area NTE is not necessarily diagnostic of NAC since there are 
competing NTE mechanisms— e.g. in-plane tilts—which will 
contribute to NTE but will drive positive area compressibil¬ 
ity under hydrostatic conditions. So in searching for layered 
NAC candidates, the strategy that emerges is to focus on sys¬ 
tems where any low-energy in-plane modes are frustrated; e.g. 
via network interpenetration, as in silver(I) tricyanomethanide 
itself (Tn . 
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